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ABSTRACT 

Fast and slow magnetosonic shock formation is presented for stationary and axisymmetric 
magnetohydrodynamical (MHD) accretion flows onto a black hole. The shocked black hole ac- 
cretion solution must pass through magnetosonic points at some locations outside and inside 
the shock location. We analyze critical conditions at the magnetosonic points and the shock 
conditions. Then, we show the restrictions on the flow parameters for strong shocks. We also 
show that a very hot shocked plasma is obtained for a very high-energy inflow with small number 
density. Such a MHD shock can appear very close to the event horizon, and can be expected as 
a source of high-energy emissions. Examples of shocked MHD accretion flows are presented in 
the Schwarzschild case. 

Subject headings: accretion, accretion disks — black hole physics — MHD — relativity — shock waves 



Introduction 



In connection with the activity of active galactic nuclei (AGNs), compact X-ray sources and inner engines 
of gamma-ray bursts (GRBs), we consider accreting plasmas onto a black hole. The observed huge energy 
outputs are mainly originated from the gravitational energy released from infalling matters. When a black 
hole is rapidly rotating in an external magnetic field, we can also expect the release of the rotational energy 
of the black hole by the electromagnetic interaction between the black hole's spin and the magnetic field (see 
Blandford & Znajek (1977) in the force-free hmit; Takahashi et al. (1990), hereafter TNTT90; Hirotani et al. 
(1992) for ideal MHD fiow). The global magnetic field around a black hole should be enhanced by infalling 
magnetized plasmas, and the magnetic field generated by the dynamo motion of the surrounding plasma 
can also interact with the spin of the black hole. The "black hole magnetosphere" system is composed of 
a rotating black hole, surrounding plasmas and magnetic fields, and considered as a central engine of AGN 
and GRB activities. It would output energy in the form of high-energy radiation (X-rays and 7-rays) and 
highly-accelerated plasma outfiows (relativistic jets and winds). In both energy output processes, the study 
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of plasma behavior in the magnetosphere is the key to various problems. This situation is very similar to 
the activity of the solar magnetosphere (i.e., the emission of X-rays and the generation of solar wind). In 
this paper, we investigate the basic plasma physics of the vicinity very close to a black hole including the 
magnetosphere. 

The magnetosphere is assumed stationary and axisymmetric, and the ideal MHD approximation is 
assumed. Then, the MHD plasmas stream along magnetic field lines in the magnetosphere, and five field 
aligned parameters exist; that is, the total energy, total angular momentum, number flux per magnetic flux 
and entropy of the MHD flow, and the angular velocity of the magnetic field line (see next section). The 
flows accelerated by gravity pass through magnetosonic points, where the poloidal component of the flow 
velocity becomes the magnetosonic wave speed. To fall into the black hole, the accreting plasma ejected 
from the plasma source with low-velocity must pass through the fast and slow magnetosonic points and the 
Alfven point, which are singular points in the basic equations. The physical flow solution along a magnetic 
fleld line is restricted by the regularity conditions at these singular points (Takahashi 2002a). That is, the 
conditions restrict the possible ranges of physical flow parameters. 

Hereafter, we will discuss the A^EHD shock in accretion onto a black hole. Rclativistic MHD shocks arc 
discussed by, e.g., Ardavan (1976) and Lichnerowicz (1967, 1976). We apply their formalism to a black hole 
magnetosphere in Kerr spacetime (for the slow magnetosonic shock, see Takahashi et al. (2002), hereafter 
TRFT02; Takahashi (2002b); Rilett (2003)). The trans-magnetosonic accretion flow, which is ejected from a 
plasma source and passes though the first magnetosonic point, results in the shock on the way to the event 
horizon, and the postshock sub-magnetosonic flow passes through the second magnetosonic point (located 
inside the shock front) again. Such a MHD shock accretion solution must satisfy the critical conditions at 
both the first and second magnetosonic points and the shock conditions. Thus, we must adjust the five 
parameters of the flow to the physically acceptable shocked accretion solution with multiple magnetosonic 
points by considering both the critical conditions at the magnetosonic points and the jump conditions at the 
shock front. 

In a plausible black hole magnetosphere, global magnetic field lines are generated by the surrounding 
plasma (e.g., an equatorial disk and its corona). We expect that (1) the magnetic field lines shaped like 
loops would be distributed around the inner-edge region of the disk, (2) the disk surface and the black hole 
are connected by the magnetic field lines, where on the disk side of the field lines the footpoints can be 
anchored to the disk surface at the location of several times the inner-edge radii and on the black hole side 
the field lines can connect the event horizon from its pole to equator, and (3) the open magnetic field lines 
connecting the outer region of the disk surface are extended to distant regions (van Putten 1999; Tomimatsu 
& Takahashi 2001; Li 2002; Uzdensky 2005). The ingoing (or outgoing) flow is ejected with non-zero velocity 
from the footpoint on the disk surface, which is a plasma source and is called the "injection point" . Then, 
along the disk-black hole connecting magnetic field lines, the MHD ingoing flows fall into the black hole, 
while the outgoing winds would stream along the open magnetic field lines and the loop-like magnetic field 
lines would enclose static plasmas. (For a high accretion rate accretion disk system, the plasma can fall into 
the black hole radially from the disk's inner-edge. In this case, the loop-like configuration of the magnetic 
field lines would disappear.) By considering the disk black hole magnetic field lines, non-equatorial accretion 
fiow that falls into the black hole from the high-latitude region is possible. If a MHD shock arises on the 
non-equatorial inflows and forms a very hot plasma region by the shock, we can expect high-energy emissions 
in the magnetosphere, which would be distinct from the emissions from the equatorial accretion disk. Such 
high energy emissions from the hot plasma would directly (or indirectly) carry informations for the strong 
gravitational field to us. 
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In our numerical demonstrations for MHD shocked ingoing flows, we should solve the realistic mag- 
netic stream function (Nitta, Takahashi, feTomimatsu 1991) for magnetic field lines connecting the disk and 
black hole. Note that at the shock front magnetic field lines bend toward toroidal and/or poloidal direc- 
tions. However, in general, it is hard to solve self-consistently the magnetic structure of the magnetosphere. 
So, for simplicity, without a realistic magnetic stream function, wc would only discuss the ingoing trans- 
magnetosonic flows on a conical magnetic stream function. We assume that at the shock front the magnetic 
field lines bend only to the toroidal direction. This assumption would be valid near the black hole (at least 
inside the inner-light surface). Near the injection point, of course, the conical magnetic fields would not 
be realistic. This is because the conical magnetic field lines do not connect to the equatorial disk surface; 
that is, we cannot define the footpoints. (Note that we can consider coronal gases distributed above the 
disk surface as a plasma source. In this case, conical magnetic field may be probable.) Furthermore, in the 
demonstrations in §4, we only treat the inflow streaming near the equatorial plane. However, we can expect 
that the qualitative picture is not drastically changed along the magnetic field lines (2) where we leave the 
equatorial region. 

In this paper, we discuss the general relativistic effects on the streaming MHD plasma in the black 
hole magnetosphere: trans-magnetosonic accretion and MHD shock formation. The accreting flow must 
be super-magnetosonic at the horizon, and the flow injected from the plasma source with low-velocity must 
pass through the magnetosonic points (the slow magnetosonic point, the Alfvcn point, and fast magnetosonic 
point). In §2, we present the general relativistic MHD flows. We introduce the field-aligned flow parameters 
and trans-magnetosonic solutions. The MHD shock formation is studied and the shocked accretion solutions 
are presented in §3. We explore various types of shocked solutions. The properties of the MHD shock is 
discussed in §4. We can obtain a very hot plasma region near the event horizon for the MHD shock of a 
high-energy inflow with small number density, where the energy conversion from the kinetic energy to the 
magnetic energy is restricted because of the black hole boundary conditions on the toroidal magnetic fleld. 
Summary and conclusions are given in §5. 



2. Critical Conditions for Trans-Magnetosonic Flows 

We consider stationary and axisymmetric ideal MHD accretion flows in Kerr geometry. In Boyer- 
Lindquist coordinates with the c = G = 1 unit, the metric of a Kerr black hole is given by 

ds^ = (l - ^) dt^ + (1) 
E ^ A 

where A = — 2mr + , T, = + cos^ 9, A= (r^ -|- a^)^ — a^A sin^ 9, and m and a denote the mass and 
angular momentum per unit mass of the black hole, respectively. The ideal MHD condition is u'^Fap = 0, 
the particle conservation law is {nu°')-a = 0, where u" is the fluid 4-velocity, n is the number density of the 
plasma, and = A,y-^ — A/^.^, is the electromagnetic tensor {A^ is a vector potential). The equation of 
motion is T^^ = 0. The energy- momentum tensor is given by 



(2) 



where fj, = {p + P)/n = rripart + (-P/^)r'/(r — 1) is the relativistic enthalpy, P is the gas pressure, p is 
the total energy density, F is the adiabatic index, mpart is the particle's mass, and F'^ = F^j^F^'^. The 
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electromagnetic field tensor Fap satisfies Maxwell's equations, F^-^.p = —A-Kj"" and *F"'^ .^p = 0, where 
*Fai3 = {l/2),/^€ai3'ysF^^ is the tensor dual to FajS and is the electric current density. The determinant 
of the metric is g, and = Ssin^. We also assume the relativistic polytropic relation P = Kp^ 

(Tooper 1965), where po = nmpa.Tt is the rest mass density, which is related to the total energy density by 
p = Pq + P/{T — 1), and ii" is a constant along the stream line for an ideal gas. The magnetic and electric 
fields seen by a distant observer, which are expressed in the Boyer-Lindquist coordinates, are defined by 
Ba = *Fafjkf^ and Ea = Fapk^, where k" = (1, 0, 0, 0) is the time-like Killing vector. 

In a stationary and axisymmetric magnetosphere, we can define magnetic field lines as \E'(r, 9) = constant 
lines, where A^) is the magnetic stream function. The plasma streams along the magnetic field line with 
five constants of motions (see Camenzind 1986a): the angular velocity of the field lines, riF(^') = —Ftr/F^r = 
—Fts/F^g, the particle number flux per unit magnetic flux, r]{'i>) = —nu^Gt/B^ = —nu^Gt/B^ = nu*{Q — 
Q.F)p'i,/ B^, where Gt = gtt + gt<p^F and Cl = u'^/u^ is the angular velocity of the plasma, the total energy of 
the magnetized flow, E{^) = fxut — ClFB^/{4:Trr]), the total angular momentum, L{^) = —p.u^ — B,p/{A-Kri), 
and the entropy >S'(5'), which is related to K. Then, we also find the relations Er = ^y—gftpB^/Gt, 
Eo = —^/^^pB^/Gt and Bt = Et = E^ = 0. From the poloidal components of the equation of motion 
with five field aligned parameters, we can derive the general relativistic Bernoulli equation (the poloidal 
equation) (Camenzind 1986b, 1989, TNTT90) and the general relativistic Grad-Shafranov equation (the 
trans-field equation) (Camenzind 1987; Nitta, Takahashi, &Tomimatsu 1991). 

The poloidal equation can be expressed by (see, e.g., TNTT90) 

I^)\l + ul)=(^)\{a-2M')f-k] , (3) 

where a = gu + igt^Qp + g^,^,^^, f = -(G^ + Gt^/lp^^M"^ - a)], k = {g^^ + 2gt4>L + gttL'^)/pl,. G,p = 
gt,p+g<p(j>^F = g,j,4,{^F — , L = L/E and u = —gt(p/g<p(p- The relativistic Alfven Mach number M is defined 
by 

Bl Bp 

where the poloidal component Up of the velocity is defined by = u^ua {A = r,6), and the poloidal 
component Bp of the magnetic field seen by a distant (lab- frame) observer is defined by Bp = —B^BaG^^. 
The toroidal component of the magnetic field i?^ = {A/T,)F0r can be reduced to B^ = —An'qEpwf, which 
is expressed in terms of the field-aligned flow's parameters and the Alfven Mach number. The denominator 
of the fimction / becomes zero when = a, where the poloidal velocity \up\ equals the relativistic Alfven 
wave speed. To obtain the physical MHD flow, which transits from sub-Alfvenic to super-Alfvenic, we must 
require that the numerator also must be zero. Then, from this regularity condition of the function /, we 
obtain L = g^^{LOA — 0,F)/GtA (TNTT90). Note that, for given L and ftp, we can find one or two radii 
satisfying the above condition for L. We will denote this radius r = r^, the "Alfven radius", while the point 
satisfying both = a and r = ta is called the "Alfven point" . The label 'A' indicates quantities at the 
Alfven radii. Note that, even if ^ a, the function / can be zero at the Alfven radius (not the Alfven 
point); the toroidal magnetic field i?^ becomes zero, while at the Alfven point B^ has a non-zero value. We 
will denote the zero toroidal field location on a magnetic field line, the "anchor point" (Punsly 2001). For a 
streaming MHD plasma, B^ changes the sign across the anchor point. 

In the case of a rotating black hole magnetosphere, there are two light surfaces given by a = 0: r = rjf (4") 
and r = rL"*(^'), where the rotational velocity of the magnetic field line becomes the light velocity as seen 
by a distant observer. The distribution of light surfaces depends on Oi?(^), a and ^{r,6). The inner light 
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surface distributes near the event horizon and encloses it. The plasma source must be located between 
the inner and outer light surfaces. For a slowly rotating black hole case, < wr < ^f, where lor is the 
angular velocity of the black hole, and for a counter rotating black hole, Of < 0, two locations of the Alfven 
radii are possible along a magnetic field line: r = and r = r^"*(^'). These radii are always located 

between the inner and outer light surfaces. Corresponding to two Alfven radii, wc can find two groups of 
the trans- Alfvenic MHD flow solutions. In one of these the poloidal velocity equals the Alfven wave speed at 
the inner Alfven radii (called the inner Alfven point), while in the other at the outer Alfven radii (called the 
outer Alfven point). ( Note that these definitions of the inner and outer Alfven points are slightly different 
from those in TNTT90.) In these cases, there is the minimum values of Lftp (= L^-Of), for which the 
Alfven points exist in the magnetosphere. There, L = L!j_(> 0) for luh < and L = !/!(< 0) for < 
satisfy the relation dL/dr^ = 0. On the other hand, for a rapidly rotating black hole case < ftp < coh, 
the number of the Alfven point is only one between the two light surfaces for an arbitrary L value. Note 
that the Alfven points are located on the light surface when L^p = 1; this corresponds to the magnetically 
dominated limit case. 

The critical conditions at the fast and slow magnetosonic points are given by the differential form of 
the poloidal equation (3), which is expressed as : 

{InUpY =J\f/V , (5) 

where 

^ = (j^y {[niM^-a)C!^ + MU^]ilnBpy 

+ 1(1 + Cl) [M\M^ - a)k' - Qa'] | , (6) 

V = {M^-af[{Cl^-ul){M^-a) + {l + ul)M^A^n-^] (7) 

with A^ = e^ + ak = /^(M^ - af, n = ae'^- 2e^M'^ - fcM^, Q = ae^ - Se^M^ - 2kM'^, and e = 1 - LCtp. 
The prime (') denotes dr + {B^ /B^)d0, which is a derivative along a stream line. The relativistic sound 
velocity is given by 

aL-(§|^) =(r-i)^, (8) 

and the sound four-velocity is given by Csw = Osw/ -\/l — a,l^- 

When the poloidal velocity equals the fast or slow magnetosonic wave speed, the numerator D becomes 

zero. The regularity of the physical solution, to be satisfied by the condition A/" = 0, is also required at 
this critical location. The fast and slow magnetosonic points have X-type (saddle-type: physical) or 0-type 
(center-type: unphysical) topology on the integral curves. The condition 2? = A/" = is expressed in terms 
of four field-aligned conserved quantities. When the values of the three conserved parameters are given 
at the plasma injection point, the value of the remaining parameter is specified by the critical conditions. 
Takahashi (2002a) discussed this problem in the Kerr geometry; for example, for given fiow parameter sets of 
L, Q.F and Ccr, the r]-rci and E-ra relations were presented, where ^cr = (asw)cr is related to the entropy and 
the index "cr" means the fast or slow magnetosonic point (to indicate the fast or slow magnetosonic point, 
we can replace "cr" by "F" or "S"). The fast magnetosonic point can be located between two Alfven radii, 
when (0 <)Llflp < Lflp < L'^''''^"''''np{< 1). Here, Z;^'^'^/™ 

are defined as L satisfying the relations both 
d[a{rA)]/drA = and a(rA) = (where the inner and outer Alfven radii coincide); L™*^ is used for a; < O 
and Z™'" is used for fi^ < 0. In this case, the hydro-like accretion solution discussed by Takahashi (2000, 
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2002a) is possible. We will call such a fast magnetosonic point the middle fast magnetosonic point. The 
hydro- like solution passes through the outer Alfven point and the middle fast magnetosonic point, and then 
falls into the black hole. However, for a rapidly rotating (or counter-rotating) plasma with \L\ > | 
the hydro-like solution is forbidden and the middle fast magnetosonic point disappears. This is due to the 
efficient centrifugal barrier on the fluid, so that no solution oi V = Af = for the fast magnetosonic point 
appears between the inner and outer Alfven points. On the other hand, the magneto-like solution is possible 
for Lj_Qp < Lflp < 1 with the condition r/ < r]™'^^. The magneto-like solution passes through the inner 
Alfven point and the inner fast magnetosonic point located between the event horizon and the inner Alfven 
point, and then falls into the black hole. In this case, the fluid part of the angular momentum transports 
to the magnetic part of the angular momentum along the ingoing flow solution. The maximum value rymax 
for the inner fast magnetosonic point exists for a hot MHD flow, while r/max — > oo in a cold limit. We can 
expect that the inner fast magnetosonic point disappears for a strong MHD shock with larger ry; such a 
MHD shocked inflow without the inner fast magnetosonic point is unphysical because the radial velocity 
of MHD inflow solutions becomes zero at the event horizon, where the number density diverges. The slow 
magnetosonic points are also obtained from the same regularity condition V = M = 0. We can find the 
inner, middle and outer slow magnetosonic points for a suitable ry- value, which arc located in the ranges of 
tl < r™ < r^, < r™'^ < r^"*, r^"* < rg"* < rg"*, respectively. For the appearance of the inner or outer 
slow magnetosonic point, there is the restriction on the ry- value (see Takahashi 2002a). 

In summary, in the case of ur < ov U,p < Q, L must have the value within the range of L'^VLp < 
Lflp < 1 to obtain two possible locations of the Alfven points in the Up-r (or M'^-r) plane. The magneto- 
like accretion solution is possible when < rf^^^ and is effective for the magnetically dominated flows. The 
hydro-like accretion solution is realized when L'^Q.p < L^p < L'^^^^'^^^p, and it becomes effective for 
the hydrodynamically dominated flows. For these two types of solutions, which also depend on E, rj and S, 
the allowable parameter ranges are distinct, so that a discontinuous transition can be expected if the flow 
parameters change in a secular timc-scalc. Note that, corresponding to two locations of the Alfven point, 
two clearly distinct types of accretion solutions exist; the magneto-like and hydro-like accretion solutions. 
On the other hand, for a rapidly rotating black hole case, < Clp < loh, only one Alfven point exists, so 
that the distinction between the magneto-like and hydro-like solution is not clear. In the next section, we 
will only treat the case of Hp > luh and mainly consider the transition from the hydro-like solution to the 
magneto-like solution through the standing shock formation. 

3. MHD Shock formation in curved spacetime 

When we consider the shock formation for accreting MHD flows onto the black hole, we must require a 
fast (or slow) miiltiple magnetosonic solution, and must apply the general relativistic MHD shock condition 
between two fast (or slow) magnetosonic points. As mentioned in the previous section, the critical conditions 
restrict the allowable ranges of physical flow parameters. In such parameter ranges, we should restrict still 
more the flow parameter ranges for multiple magnetosonic shocked accretion solutions. 

3.1. The jump conditions 

Here, we discuss the jump condition for MHD flows in Kerr geometry. The particle number conservation, 
the energy momentum conservation and the magnetic flux conservation across the shock in a relativistic MHD 
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flow are (Lichnerowicz 1967, 1976; Appl & Camenzind 1988) 

[nu"]l £^ = 0, (9) 
[T"^]2 = , (10) 
[*F-^]l 4=0, (11) 

where i" is the unit vector normal to the shock front, which has only poloidal components for stationary 

and axisymmetric flows {£"^£a = —1)- The symbol [Z]l = Z\ — means the change of a certain quantity 
Z across the shock located at the radius r = rsh(^'), where the indeces "1" and "2" mean the properties of 
the preshock and postshock flows just on the shock front, respectively. From the conditions (9) and (11), 
the number flux across the shock front U = nu"£a and the normal component of the magnetic fleld (seen 
by a distant observer) to the shock front B± = B°'£a remain unchanged; that is, [U]l = and [B±]2 = 0. 
We can also obtain [E^^]^ = 0, ['^f]2 — and \rf^ — 0, where = eABE"^£^ is the tangential component 
of the electric field E". The quantity rj can be expressed as 77 = —UGt/B±. From the condition (10), the 
following vector remains unchanged across the shock: 

w'.,.<.u-Pr-^^ (§-)[§-), (12) 

where P = P+{B'^/8Tr) and ^ {F'^/2) = aB^ + Bl/pl . Prom the product W"£a, we obtain the relation 

[(M/n)C/2 +P]l = 0, (13) 

which can be reduced to 

(M/n)i - (M/n)2 - M| ^ ' 

Prom the conditions [W*]^ = and [W^''^]2 = 0, we also obtain that the constants of motion E and L are 

continuous across the shock, while entropy S, which is the fifth constant of motion, is discontinuous across 
the shock (see TRFT02). Of course, for the physically realistic shock solution the entropy must increase. 
Then, we will introduce the entropy related mass fiow rate per magnetic fiux tube defined by Chakrabarti 
(1990) 

M = K^'n.-n , (15) 

where N = l/(r — 1) is the polytropic index. The quantity Xd conserves for a shock-free flow, but it must 
increase at the shock due to the entropy generation. By using the definitions of the Alfvcn Mach number 
and the relativistic enthalpy, we can express the entropy-related accretion rate as a function of r and 
with the conserved quantities 



M 



1 / -^-1^ 
1 f 1^ 1 



1 + N 



(16) 



where from the poloidal equation the relativistic enthalpy jj can be expressed as 



Thus, we can plot trans-fast MHD solutions as A4{r,M'^) = constant (> 0) curves on the r-M"^ plane. The 
physically acceptable MHD shock must satisfy the condition < Mi < M2 ; M. = 0fora cold fiow (/U = /Zc). 
The flow in the A4{r, M^) < region is forbidden as a physical solution. 
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The location of the standing shock front and shock properties depend on the field-aligned flow param- 
eters. To discuss the shock properties, we will introduce the following dimensionless parameters (Appl & 
Camenzind 1988, TRFT02): 



? - a = > (18) 



n2^4 _ M^je-hMj) 
^ ~ niu\ Miie-hM^f' ^ ^ 

where h = g**{l — Llo). The parameter q is the ratio of the toroidal magnetic field before to that after the 
shock, and ^ is the compression ratio; these parameters are quantities seen by a distant observer, and include 
the gravitational red-shift factor and Lorentz factor for the plasma motion. We also define the plasma frame 
compression ratio A by 

Furthermore, wc arc interested in the jump of the magnetization rate and temperature of the plasma. The 
magnetization parameter a, which is defined as the ratio of the Poynting flux to the total mass-energy flux 
seen by a zero angular momentum observer (ZAMO), can be expressed as (TRFT02) 

_ B^G^ _ e- ah 

The temperature parameter is defined by 

e = fcBT/mpart = Kp'J^ = (m/mc - 1)/(1 + N) , (22) 
where T is the temperature of the fluid and fee is the Boltzmann constant. 



3.2. Shocked accretion solutions 

Now, let us consider MHD accretion solutions with a fast/slow magnetosonic shock on the r-M^ plane. 
First, we need to find two solutions of the trans-magnetosonic flow (i.e., trans-fast or trans-slow magnetosonic 
flows corresponding to the fast or slow magnetosonic shock) having the same values of E, L, rj and ^Ip- By 
the shock formation, wc can connect these two trans-magnetosonic solutions as a miiltiplc magnetosonic 
accretion solution. Then, the shocked accretion flow passes through the fast or slow magnetosonic points 
twice. Along each branch of the preshock and postshock solutions, the value of M. (or the entropy S) is 
different, A^prc 7^ A^postj where the indexes "pre" and "post" indicate the preshock and postshock trans- 
magnetosonic accretion solution. The values of Mp^e and Alpost are specified by the critical condition at each 
fast or slow magnetosonic point for the preshock and postshock solutions; that is, TWpre = -M^rcri, M^ii) 
and jVlpost = .A4(rcr2, -^cr2)) where the labels 'crl' and 'cr2' indicate the first and second fast or slow 
magnetosonic points, respectively. 

When wc require a MHD flow solution passing through the slow magnetosonic point S, the Alfvcn 
point A and the fast magnetosonic point F (hereafter "SAF-solution" ), from the critical conditions at the 
fast and slow magnetosonic point, a relation between E, L, rj and is also specified. So, a change of 
the value of one of the parameters requires the change of the values of the remaining parameters. In the 
next section, to demonstrate fast magnetosonic shocked accretion solutions, we will treat the hydro-like 
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Fig. 1. — Schematic pictures of physically acceptable accretion solutions with the fast magnetosonic shock 
(thick curves and arrows). The upstream SAF-solutions are (a) hydro- like and (b) magneto- like, and jump 
to the downstream trans-fast magnetosonic solution by the shock formation. The large circle in the Left 
panel (a) represents the plasma source, where the plasma is injected with a finite velocity. In the Right 
panel (b), the downstream solution docs not pass though the inner Alfvcn point, but connect to the event 
horizon with a finite Mach number. The vertical dotted lines are the Alfven radii. The boundaries of the 
hatched regions are A4 = curves and are labeled by "0" . The A4 value of the downstream curve with the 
label "2" is grater than that of the upstream with the label "1"; that is, the entropy of the flow is increased 
by the shock. The other partially plotted thin curves that do not connect the plasma source are unphysical 
accretion solutions. 

SAF-solution as an upstream flow solution and the magneto-like solution as a downstream flow solution for 
the fast magnetosonic shock formation. When we demonstrate a slow magnetosonic shock formation, we 
will consider that the downstream flow solution is the magneto-like SAF-solution, while the upstream flow 
solution is trans-slow magnetosonic solution. Then, we can restrict the values of Alprc and A^post for the 
MHD flow solution passing through the first and second X-type fast magnetosonic (or slow magnetosonic) 
points, respectively. To find such a SAF-magnetosonic accretion solution, for a given parameter set of ^p, 
E and L we can search for the acceptable value of by tuning the regularity conditions (I? = A/" = 0) at 
both fast and slow magnetosonic points. If a suitable parameter set is obtained, we will find at least two 
X-type fast/slow magnetosonic points on the r-M^ plane. 

Next, by applying the relativistic jump condition, we can find one (or two) shock location to give 
7VI2 = Alpost with the condition M.i = Alprc- We will demonstrate an accretion solution with the MHD 
shock. The accretion solution with a shock is obtained on the r-M^ plane by plotting A4{r, M^) = constant 
curves of the two trans-fast magnetosonic (or trans-slow magnetosonic) flow solutions and a vertical line 
connecting these solution curves at the shock location. 

Although the numbers (and their locations) of the fast and slow magnetosonic points depend on the 
flow parameters and the spin of the black hole (Takahashi 2002a), we can find the following cases of the 
shocked MHD accretion solutions: two classes of fast magnetosonic shock accretion solutions. 



FS-i : inj 




pmid 



< shock > 



F 



H 
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(a) (b) 




Fig. 2. — Schematic pictures of unrealistic accretion solutions. The upstream solutions are SAF-solutions, 

but are broken before the event horizon. In this case, the entropy (or A4) of the upstream flow labeled by 
"2" is grater than that of the downstream labeled by "1", so that the shock formation is forbidden (thick 
broken lines mean the disconnection of the upstream and downstream solutions). If the downstream curve 
is SAF-solution and connects to the plasma source and the event horizon, it is a physical accretion solution 
without a shock. 

FS-ii : inj ^ S-'d _ A- ^ f 'Zter ^ < shock > ^ FP„,, ^ H , 
and one class of slow magnetosonic shock accretion solution, 
SS : inj ^ S°"* ^ < shock > S"^'*^ ^ A'" ^ ^ H . 

In Figure la (case FS-i), the upstream flow is a hydro-like solution and the downstream postshock flow is 

a magneto-like solution. The shock formation is possible between these two trans-magnetosonic solutions, 
where the requirement M.i < M.2 is satisfied. For the physically acceptable accretion solution, the postshock 
trans-magnetosonic solution must connect to the event horizon with nonzero four-velocity, = finite. The 
upstream solution passes through the outer Alfven point A°"* and the first X-type middle fast magnetosonic 
point F'"''^ After the fast magnetosonic shock, the downstream solution passes throiigh the second X-type 
inner fast magnetosonic point F'", and then falls into the black hole H. Figure lb (case FS-ii) shows that two 
X-type inner fast magnetosonic points (inner F™ and outer F'") locate inside the inner Alfven point. The 
upstream magneto-like solution passes through the first X-type inner fast magnetosonic point FJJJjj^j., after 
passing through the inner Alfven point A'". The upstream solution can transit to the downstream trans- fast 
magnetosonic accretion solution, which passes through the second X-type inner fast magnetosonic point 
Fj^ner and connect to the event horizon H. On the other hand, in the case of (SS), the upstream solution of 
the slow magnetosonic shock connects to the outer Alfven point, but before the Alfven point the upstream 
solution can jump to the downstream SAF-solution. This downstream solution passes though the middle 
slow magnetosonic point S™"^, the inner Alfven point A™^'^ and the inner fast magnetosonic point F'", and 
then falls into the black hole. 

One may consider that an upstream accretion solution does not need to connect to the event horizon H. 
Such a solution is of course unphysical as an accretion solution, but by the shock formation this upstream 
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solution may transit to a downstream physical trans-fast MHD accretion solution that connects to the 
horizon with a finite (non-zero) four- velocity. Although the upstream SAF-solution has M. = Alpre and the 
downstream trans- fast MHD solution has M = Alpost) the M. = TVlpre curve makes a loop or a double- loop 
on the r-M ^ plane by connecting to the inner Alfven point (and then it is disconnected to the event horizon) 
as shown in Figs. 2a and 2b, while the Ai = Alpost curves enclose this loop and one of them connects to H. 
However, by considering the entropy (or M) distribution, we conclude that the entropy of the disconnected 
upstream solution is greater than that of the downstream trans- fast MHD solution; that is, Mp^e > Alpost- 
So the transition from the upstream to downstream solution is forbidden. Thus, for the fast magnetosonic 
shock formation, both upstream and downstream trans-fast accretion solutions need to connect to the event 
horizon with non-zero velocity. We should note that even if two trans-magnetosonic accretion solutions exist 
for a set of values of the field-aligned quantities, no MHD shock may generate under such values. That is, 
we may find a MHD shock of Ali > A42, but such a solution is unphysical. 



4. Numerical Results and Discussions 

Now we show class (FS-i) fast magnetosonic shock solutions, where a cold {& <C 1) preshock flow is 
assumed. The solutions are given for the equatorial plane {9 — n/2) in the Schwarzschild spacetime (a = 0). 
For simplicity, we set the shock normal to the downstream flow; that is, = {0,£^ ,0,0). Then, we obtain 
the relations Bl = B2 = B± and Ef = E2 = E\\ . Furthermore, without the trans-field equation the poloidal 
magnetic field is assumed to be the radial configuration denoted by B^ = BoGt/^ and = ( that is. 
Bp = Bp/ Bo = l/\/AS) between the plasma source and the event horizon, where So = constant (7^ 0) 
and the value should be determined at the plasma injection point. Some special radii for the MHD flow are 
given by the field aligned parameters. The locations of the light surfaces are given by the Clp value, and 
the locations of the Alfven radii are given by the ^lF and L values. In addition, setting F = 4/3, the values 
of E = E/fjbc and 7) = iwUcri/BQ with the values of Of and L determine the locations of the fast and slow 
magnetosonic points. Then, solution curves M"^ = M'^{r;Ai) are obtained on the r-M"^ diagram, and the 
final parameter M is used to specify an acceptable trans-magnetosonic solution for accretion. ( The solution 
curves also include an outgoing physically acceptable solution that can take a different value of M. with the 
accretion solution. The most remaining solution curves are unphysical. ) 

By assuming a cold preshock flow {M.i = M.pie = 0), the shock condition (13) can be reduced to 



I pre 

2 /.\2 ..,_p 



^\ M',+lE'f, = i^ \ Mi + {f,B,M,f-' 



M2 



+ \E^fi . (23) 



This equation refers to the jump of the Mach number across the shock front. When we obtain the two 
trans- fast (or trans-slow) magnetosonic solutions = Mp^g(r) and = M'^^^{r) by tuning both E and 
fj values with the requirement oi AAi = 0, we can apply the solutions to equation (23); that is, we can 
set Ml = Mpjg(rsh) and M| = MpQg^(rsh)- Then, we obtain the location r = rgh satisfying the condition 
M.2 — M.post = 0. Note that for an acceptable shocked accretion solution we must require the condition 
M.'i > M.I, while the other field aligned parameters Of, E, L and 77 are conserved across the shock. 

Figure 3a shows the square of the Alfven Mach number of the shocked accretion solution vs the radial 
distance. The solution, which passes through the fast magnetosonic points twice, is represented as two M. = 
constant curves and a vertical arrow (the fast magnetosonic shock). The upstream SAF-curve follows the 
hydro-like solution, although the slow magnetosonic point is absent on this diagram because of the cold 
approximation. Then, at the shock location, it transfers to the downstream curve which is the magneto-like 
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Fig. 3. — The trans-fast magnetosonic accretion solution with a fast magnetosonic shock (thick black curves); 
(a) The square of the Alfven Mach number A'P vs. radius r/m and (b) the toroidal component of the magnetic 
field Bfjj (= B^/Bq) vs. radius r/m relation are shown with BqM = constant. After passing through the 
outer- Alfven point (A) and the middle- fast magnetosonic point (F™"*), the ingoing flow makes the fast 
magnetosonic shock near the event horizon. The postshock flow passes through the inner-fast magnetosonic 
point (F'"), and then falls into the black hole (r = 2m). The gray curves also show B^Xi — constant, 
but these are unphysical solutions. The flow parameters are given by ?7 = 0.0108, E = 3.5, L/m = 4.1, 
m^F = 0.14503, r = 4/3, a = and 6* = n/2. We obtain B0M2 = 0.05134 for the postshock flow, while we 
give i?o-^ 1 = for the preshock flow. 

solution, and it becomes sub-fast magnetosonic. The shocked flow passes through the fast magnetosonic point 
again to fall into the black hole. Figure 3b shows the toroidal component of the magnetic field {B^ = B^/Bq) 
for the shocked solution as a function of radial distance. As the plasma falls inward, the trailed-shape of 
the magnetic field line {B^ < 0) changes to the leading-shaped one {B^ > 0) at the Alfven radius (the 
anchor point; not the Alfven point), where = 0. The strength of the postshock toroidal magnetic field 
\Bij)\ increases across the fast magnetosonic shock. In Figure 3 the preshock flow solution is plotted as a 
bold curve starting from the outer light surface with zero-velocity, where the plasma rotates toward the 
(/(-direction with the speed of light. So, a realistic flow should start from a location somewhat inside the 
outer light surface. In the case of (FS-i), the preshock solution should connect to the plasma source at some 
location between the outer-Alfven point and the outer light surface. In Figure 3, gray curves also show A4 = 
constant curves, but they do not pass though the magnetosonic points, or do not connect to the plasma 
source, except for the upstream = curve of r/f < r < 7-sh- Note that the upstream Ai ~ Q curve 
directly connects the plasma source to the event horizon, where the Mach number has a finite value at the 
event horizon. This curve is also a physically acceptable SAF-solution without the MHD shock. Here, we 
cannot discuss which of the two accretion solutions; the shocked inflow or the shock-free inflow, is selected 
as a realistic accretion flow. To answer this question, the stability analysis for shocked MHD flows would be 
necessary, but that is beyond the scope of our current paper. In the following, we will go ahead to further 
explore the shocked MHD accretion flows. 
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Fig. 4. — The relations between the shock location rsh/m and (a) the energy E, and (b) particle number 
flux fj for MHD shocked accretion flows onto a black hole (FS-i solutions). The flow parameters are given 
by L/m = 4.1 (thick solid curves) and 4.15 (thick broken curves), mflp = 0.14503, a — 0, and 9 = n/2. 
The branch inside the inner Alfven radius marked by "□" shows the fast magnetosonic shock, while the 
branch between the radius marked by which is located inside the middle-fast magnetosonic point, and 
the inner-Alfven radius shows the intermediate shock. The curves labeled by "Fp" represent the radii of 
the inner/middle fast magnetosonic points. The lines labeled by "Ap" represent the radius of the Alfven 
point, and the lines with "Ar" indicate the Alfven radius. The solid and broken lines/curves correspond to 
L/m = 4.1 and 4.15, respectively. 

The parameter search for the shocked flow solutions is a standard approach to correctly understand 
the shock formation in a curved spacetime (TRFT02) and to estimate the activities of a black hole engine. 
Figure 4a shows the relation between the shock location r^h and the total energy E for the class (FS-i) 
solutions, which are obtained by solving equation (23) with M2 — Alpost = 0. Figure 4b shows the relation 
between the shock location Vgh and the particle number flux per the magnetic flux ry. The inverse of this 
quantity exhibits almost the same behavior as E. For a physically acceptable shocked MHD accretion 
solution onto a black hole, the total energy E has the minimum value £'min under a certain parameter set 
of Hp, L and a. It seems that the energy E diverge at r = r^™ (where rjE has a finite value as seen in 
Fig. 5), but we stop the calculation &i E — 100 (marked by o). Note that E 00 (and — > 0) with 
Lftp 1 corresponds to a force-free limit, where the Alfven radius shifts to the light surface; but in this 
demonstration Lilp = 0.6; that is, it is not force-free. The possible shock location for physically acceptable 
shocked accretion flows has the minimum radius J"™" ^^'^ maximum radius r"^^^ (marked by •), where 
these radius are located between the middle- and inner-fast magnetosonic points; that is, rF2 < '"sh'" ^'^^ 
r'^^ < rpi- In Figure 4, the radii of the inner-fast and middle-fast magnetosonic points weakly depend on 
the value of E or 77. The inner Alfven radius and the outer Alfven point are plotted by the vertical lines. 

Figure 5 shows the value of the total energy flux per magnetic flux tube rjE at a fast magnetosonic 
shock location. At a shock location, the value of rjE for a shocked flow increases with increasing L, while the 
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Fig. 5. — The relations between the shock location rgh/m and the energy flux fiE for a (FS-ii) solution. The 
multiple shock locations are shown in the magnified Right panel (b). The chosen parameter sets are the 
same as in Fig. 4. 




Fig. 6. — The relations between the shock location rsh/m and the entropy generation -BoA^2 by the shock. 
The chosen parameter sets are the same as in Fig. 4. 
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value of E increases and the value of rj decreases (see Fig. 4). Then, the fast magnetosonic shock with larger 
energy flux is obtained for the inflow with smaller E (or larger rj) and larger L. Figure 6 shows that the 
relation between the shock location and the entropy related mass flow rate per magnetic flux tube M.2 for the 
postshock flows. For the shock located near the inner-fast magnetosonic point the lager entropy generation 
is obtained, while the energy flux becomes small at there. With these flow parameter sets, we may expect a 
stronger shock formation and/or larger energy release from the hot plasma generated by the shock formation. 
We will discuss this problem later. Note that a maximum value of L exists for the hydro- like solution; The 
hydro-like solution is necessary for the middle-fast magnetosonic point to achieve the class (FS-i) solution 
(Takahashi 2002a). 

When a certain value of the total energy E (or rjE) is given, one or two shock locations are obtained 
(see Fig. 5b to confirm the multiple-shock locations). When we choose another value of E (or r]E), different 
values of t] and M.2 are determined from the critical conditions at the magnetosonic points. The number 
of physically acceptable shock location depends on the total energy E (or rjE). For example, the number is 
one for E^, < E < -Bmax, and two for i?niin < E < E^, (a multiple-shock), where E^, is given at the outermost 
shock location r^'^^. When we find the multiple-shock locations for a given parameter set, we should discuss 
the stability of MHD shocks at each shock location to answer the question of which shock location should be 
selected as an accretion solution. However, that task is beyond the scope of the present paper. An important 
point is that the maximum and minimum values (labeled by 'max' and 'min') exist; that is, the condition 
for the shock formation is limited by some ranges of the parameter values. 

If a shock generates between the inner Alfven radius and the possible outermost shock radius, which 
is a segment divided by "□" and "•" in Figs. 4-9, the postshock flow becomes sub-Alfvenic. Such a 
shock is called the "intermediate" magnetosonic shock (the stability of intermediate shocks is discussed by 
Hada (1994)). In the case of an intermediate shock, the post shocked flow must pass through both inner 
Alfven point and inner- fast magnetosonic point. In Figures 4 and 5, we see that the multiple-shock locations 
are caused in the intermediate shock region. [ For a hot MHD inflow streaming along an area-expanding 
(diverging) magnetic fleld line of Bp oc (l/\/AS)r~^, which is converging along ingoing flows, the possible 
multiple-shock locations extends to the fast magnetosonic shock region (Goto 2003). ] For smaller values of 
E (or larger rjE), multiple-intermediate shock solutions are possible, while an accretion flow with a larger 
value of E (or smaller r]E) produces a fast magnetosonic shock. 

Figure 7 shows the plasma frame compression ratio to be A < 1.7 for L/m = 4.1 and A < 1.8 for 
L/m = 4.15; the compression ratio ^ seen by a distant observer is somewhat small because the factor 
wl/wi < 1 weakens the shock strength. Note that the Lorentz factor m* includes the effects of strong gravity 
(gravitational red-shift) and the Lorentz boost by the rapid motion of the plasma around the black hole. 
The strongest shock (the largest A value) is generated around the inner Alfven radius, where the difference 
in the Alfven Mach number between the prcshock and postshock trans-fast magnetosonic solutions is also 
the largest. Although the number density ni of the cold preshock flow is inversely proportional to Mf, the 
function m = ni(r) has the minimum value around the middle-fast magnetosonic point. 

When the shock is generated just on the inner Alfven radius that is the anchor point for the upstream 
flow, we see the switch-on shock with \l/q\ = (see Fig. 8a), where (-B^)i = for the preshock flow and 
{Bcf,)2 7^ for the postshock flow. In the fast magnetosonic shock of the class (FS-i), the inner Alfven radius 
is not the Alfven point for the upstream flow solution but it is the Alfven point for the downstream flow 
solution. The upstream flow becomes hydro dynamical {E = fxut, L = —fm^) just on the anchor point. For 
the positive q, as the ratio of the preshock to postshock toroidal magnetic fields 1/q increases in magnitude, 
the magnetization parameter of the postshock flow also increases (see Fig. 8b). The property g > 1 is a 
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Fig. 7. — The compression ratio A vs. shock radius relation is shown. The chosen parameter sets are the 
same as in Fig. 4. 




Fig. 8. — (a) The inverse of the toroidal magnetic field ratio across the shock 1/qvs. shock radius relation is 
shown, while (b) the magnetization parameter a for the pre- and post-shock inflow vs. shock radius relation 
is shown. The chosen parameter sets are the same as in Fig. 4. 
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Fig. 9. — (a) Shows the temperature parameter Q vs. shock radius realtion, while (b) shows the jump of the 
radial four-velocity Au^ vs. shock radius relation. The chosen parameter sets are the same as in Fig. 4. 

characteristic of the fast magnetosonic shock. The negative q corresponds to the intermediate MHD shock, 
where the sign of the toroidal component of the magnetic field is reversed. In the class (FS-i) solution, 
the Poynting flux directs inward across the fast magnetosonic shock; that is, the electromagnetic energy is 
carried to the black hole by the preceding rotation of the anchor point. However, for the intermediate MHD 
shock, the Poynting fluxes direct outward (inward) for the preshock (postshock) flows; in the intermediate 
MHD shock the magnetic fleld line is flipped over at the shock normal. 

As seen in Fig. 3, when the shock is approaching the inner-fast magnetosonic point, the gap between 
Ml and M| becomes small. One may expect that only some fraction of kinetic energy of the preshock 
flow converts to thermal energy of the postshock flow. In fact, the strength of the shock (the compression 
ratio A) is weakened; this means that the ratio of the preshock and postshock 4- velocities (M^/u2)sh becomes 
small. To confirm this, in Figure 9a, the temperature of the postshock flow is shown as a function of the 
shock radius. Figure 9a (see also Fig. 6), however, indicates that the plasma becomes hotter for the fast 
magnetosonic shock generated near the inner- fast magnetosonic point. This is because for the shock located 
near the inner-fast magnetosonic point the plasma's energy per rest-mass energy E is very large and the 
number flux per magnetic tube ry is small. Furthermore, as seen in Figure 9b that shows the jump of the 
radial 4- velocity Au"^ = (Mi)sh — (ui)sh at the shock location, the gap Au*" becomes large near the inner- fast 
magnetosonic point even if the ratio {ui/u2)sh is small (where the value of A has a small value), and a 
considerable amount of plasma kinetic energy converts to the thermal energy. The plasma kinetic energy, 
of course, is converted to the magnetic energy by increasing i?^. However, for the fast magnetosonic shock 
generated near the event horizon, the strength of {B^)2 is ^ [oc {luh — ilF)B^], which is determined 
by the boundary condition at the event horizon (see Fig. 3b) and is independent of itj^ and E (and r}). 
We can say that the toroidal magnetic field B^ of the stationary shocked accretion flow is restricted by the 
general relativistic effect (i.e., the hole's boundary condition), although the information of the horizon cannot 
propagate to the shock front. Then, the magnetic energy generation cannot dominate in the shock process. 
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Fig. 10. — A MHD accretion flow with a slow magnetosonic shock is shown: (a) The square of the Alfven 
Mach number and (b) the toroidal magnetic field vs. radius r/m (thick black curves). The preshock curve 
does not connected to the horizon, but it is physically acceptable through the shock. The flow parameters 
are given by ^ = 0.00521, E = 7.012, L/m = 4.5, mVLF = 0.14303, V = 4/3, a and 6* = ■n/2. Wc obtain 
— 3.591 for the postshock flow, while we give BqMi = for the preshock flow. 

Thus, even if the shock is weak, the high temperature postshock inflow is obtained. We can also see that 
the temperature of the postshock inflow is rising as the shock location approaches the event horizon. So, a 
cone-like hot plasma column would be built on the event horizon. The high-energy emissions are expected 
from this hot plasma column with rn < r < Tgh- 

Figure 10 shows an accreting flow with the slow magnetosonic shock. Here, the upstream flow solution is 
cold, so that the injected plasma from a plasma source is super-slow magnetosonic. (Note that the upstream 
super-slow magnetosonic solution does not necessarily have to pass through the middle-fast magnetosonic 
point.) After the slow magnetosonic shock, the ingoing flow becomes hot, and passes through the slow 
magnetosonic point, the Alfven point and the fast magnetosonic point in this order, and then falls into a 
black hole. Although, here the slow magnetosonic solution is generated near the outer Alfven radius, the 
jump of the Mach number of the slow magnetosonic shock is larger than that of the fast magnetosonic shock 
(see Fig. 3); that is, the compression ratio becomes higher. The preshock toroidal component of the magnetic 
field decreases with decreasing radius, and at the slow magnetosonic shock it decreases catastrophically. After 
the slow magnetosonic shock, the trailed-shape of the magnetic field line changes to the leading one. The 
detailed properties of the slow magnetosonic shocks are discussed by TRFT02 and Takahashi et al. (in 
preparation) . 
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5. Summary and Conclusion 

We have formulated the MHD shock conditions in Kerr geometry and have discussed MHD accretion 
flows onto a non-rotating black hole with magnctosonic; shock formation. We find that a very hot plasma 
region can form close to the black hole. To realize these flows, solutions with multiple magnctosonic points 
are necessary, where the critical conditions at the magnctosonic points are required for the acceptable set 
of the field-aligned parameters. Furthermore, the MHD shock must be located somewhere between the 
two magnctosonic points, where the shocked plasma must satisfy the jump conditions. Although the flow 
parameters should be specified at a plasma injection source (e.g., the surface of an accretion disk or its 
corona, etc), these parameters for the acceptable MHD accretion flows are specified at some points located 
near the event horizon by the above conditions. Then, the shocked accretion phenomena would give us 
information about the plasma sources and the magnetosphere (the inflowing region) in a curved spacetime. 

A strong MHD shock with large fluid compression is obtained when the Mach number gap between 
the preshock and postshock flows is large; such a shock generates around the inner- Alfven radius. We find 
that this situation is realized through the transition from the hydro-like preshock flow to the magneto-like 
postshock flow (type FS-i shocked MHD accretion flow). In the case of a rapidly rotating black hole with 
< < ^Hi this situation with the hydro- like and magneto- like solutions is lost. However, this is possible 
for the type (FS-ii) shocked MHD accretion flows. In that case, one may not be able to expect strong MHD 
shock formation with a large compression ratio because the Mach number jump is not so large as in the 
(FS-i) case. However, a very hot plasma that would be caused by the entropy generation at the shock with 
^sh ^ would be expected as in the case of (FS-i) shocked accretion flows. 

The Lorentz factor of relativistic jets observed in AGNs is ^ 10 - 100. Therefore, we expect, as an 
origin of the jet, that the ejected plasma wind from the disk surface has the total energy of ~ 10 - 100. 
The ejected high-energy plasma streams along the magnetic field lines in the magnetosphere, and forms both 
outgoing winds and ingoing winds, where the separatrix surface, which separates infiow and outflow regions, 
depends on the magnetic field distribution in the magnetosphere. Although the outgoing wind outside the 
separatrix surface is accelerated and would make a relativistic jet, some part of the ejected wind inside the 
separatrix surface streams toward the black hole because of its strong gravity. Such a wind can produce 
a MHD shock discussed in this paper. The shocked plasma with i? ~ 10 - 100 will then make a very 
hot plasma region. Thus, we can expect that a very hot plasma region is located near the event horizon 
as a source of high-energy radiation (X- and 7-ray emissions). Of course, most radiation fluxes would fall 
into the black hole by the gravitational lens eff'ect, and the energy of the outgoing radiation is lost by the 
gravitational red-shift effect. However, a huge energy radiation generated by the shock would modify the 
distributions of the plasma and magnetic field geometry around the black hole. For example, some kind 
of plasma instabilities or dynamical phenomena, which will include the information on the strong gravity, 
may be caused by the strong radiation. Such plasma phenomena taking place between the shock and the 
inner-fast magnctosonic point can propagate outward, and may influence the upstrem flow and further the 
plasma source. In a magnetized disk-black hole system, "hyperaccretion" onto a slowly rotating black hole 
is proposed as a model for short gamma-ray bursts (van Putten and Ostriker 2001). If a very hot plasma 
region by the MHD shock is generated along the hyperaccretion flow, a fraction of accreting plasma may be 
blown away as a GRB jet by some plasma instabilities or dynamical phenomena. 

When we are interested in a dense accretion flow (with a larger 77 value) that leaves the disk surface 
without large initial velocity, the expected energy would he E < 1. Although we have explored MHD shocks 
for larger E, we also tried to search for the minimum energy of shocked MHD accretion flows (FS-i) under 
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wide ranges of flow parameters (including the variations of L, Q,f and 9 values). However, we did not flnd a 
MHD shock solution for £^ < 1 in our limited parameter search. More systematic analysis will be presented 
in our subsequent paper (Fukumura et al. in preparation). 

In the case of quiet black hole accretion system {E < 1), no MHD shock may be obtained near the 
equator. However, the MHD shock can be generated along a disk-black hole magnetic field line, and the hot 
plasma region can form in the high-latitude region of the black hole. Then, we expect that as a signature of 
AGN activities X- and 7-ray emissions from the off-equatorial shocked hot plasma are observed dircx'tly to 
us and would be capable of locally illuminating the underlying accretion disk. This could photoionize iron 
atoms in the disk, causing subsequent iron fluorescence observed in many Seyfert nuclei and Galactic BH 
candidates. Our MHD shock model thus can be a probable local X-ray radiation source in these systems. 
Although the estimate of the emergent X-ray spectrum (at the shock) is important, it is beyond the scope 
of our current investigation. 

We would like to thank Akira Tomimatsu for his helpful comments. This work was supported in 
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A. Shock Condition in Plasma Comoving Frame 

In this paper, we use the expression with the lab-frame magnetic fleld to analysis the MHD shock 
conditions. Here, we show the relations between the lab-frame magnetic field and the plasma comoving 
frame magnetic field. 

The magnetic and electric fields in the plasma comoving frame arc defined by = * Fa^u" and = 
FaffU^ , which satisfy Cq, = for the ideal MHD plasma. Then, the homogeneous Maxwell equation can 
become {u°'b^ — u^b°')-a = 0, which means the magnetic flux conservation, and equation (2) becomes 

T^"' = wu"u^ - Pg"'^ - —b°'b'^ , (Al) 

where w = nfi+ib'^/An), P = P+ib'^/Sn) and 6^ = -b'^ba- The lab-field denoted in the Boyer-Lindquist 
coordinates are related to the comoving magnetic field 6" by 

Br = V^F'f'^ = Gtbr/e , (A2) 
Be = V^F^'I' = Gtbe/e , (A3) 
B^ = ^F^"" = {Gtb^-G^bt)/e , (A4) 

where e = {E—ClpL) / fi. Prom the shock conditions, the following scalar U and two vectors remain unchanged 
across the shock (Lichnerowicz 1967) : 

U = nu± , (A5) 
V" = b±u" - u±b" , (A6) 
W" = u±wu" - Pr - {l/ATr)b±b" , (A7) 

where u_\_ = u"£a and b_\_ = b^ta (the index "_L" indicates the normal component of the vector to the shock 
front). Note that V"£a = 0. Furthermore, from the product V"Wce we obtain the relation [/u6j_]2 = 0, which 
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is equivalent to [B±]2 = 0, where B± = B^ia- Thus, the qantity B± ( or iJ,b±) in a curved-space is a shock 
invariant. Note that the normal component of the comoving-frame magnetic field changes across the shock, 
while in the non-relativistic limit we find n = Hc = Wpart, so that b± remains unchanged. 
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